Levy Statistics for Random Single— Molecule Line Shapes in a Glass 
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We demonstrate that the statistical behavior of random Une shapes of single tetra-tert- 
butylterrylene chromophores embedded in an amorphous polyisobutylene matrix at T = 2 K, is 
described by Levy statistics as predicted theoretically by Barkai, Silbey, and Zumofen [Phys. Rev. 
Lett. 84 5339 (2000)]. This behavior is a manifestation of the long range interaction between two 
level systems in the glass and the single molecule. A universal amplitude ratio is investigated which 
shows that the standard tunneling model assumptions are compatible with the experimental data. 
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Experimental advances Q have made it possible to 
measure the fluorescence of a single molecule (SM) em- 
bedded in a glass. Because each individual molecule is 
in a unique static and dynamic environment, the fluores- 
cence of chemically identical SMs varies from molecule 
to molecule. In this way the molecules serve as lo- 
cal reporters on the dynamics and statics of the host 
glass. SM experiments have beenperformed both in low- 
temperature glasses 0; S S IE and recently close to 
the glass transition temperature |3 ■ For low-temperature 
glasses the fundamental question is: Is the standard tun- 
neling model valid for glasses? Related questions are how 
to analyze the complex line shape behaviors of SMs in 
glasses, and what do their random behaviors teach us on 
the SM-glass system. 

The standard tunneling model ^ was suggested in the 
early seventies to explain universal features of glasses; 
for example many glasses show a heat capacity which 
is nearly linear in temperature. At the center of this 
phenomenological model is the concept of the two-level 
system (TLS). It is assumed that at cryogenic temper- 
atures excitations in glassy materials are two-level tun- 
neling systems whose energies and tunneling matrix el- 
ements are randomly distributed. More recently, Geva 
and Skinner modeled behaviors of SMs in glasses 
based on the standard tunneling model. Orrit and co- 
workers 0, |2] used the fluorescence of single terrylenc 
molecules in the polymer polyethylene to obtain the flrst 
direct experimental proof that two level systems actually 
exist in an amorphous material. More recent experiments 
PI revealed behaviors not compatible with the standard 
model for 21 out of 70 molecules; e.g., a molecule coupled 
to a three level system ^ . We note that a fundamental 
first- principle understanding of TLSs is still missing |llj| , 
although numerical simulations |l2j| give some evidence 
on the microscopic nature of a few of these entities. 

There are many open questions concerning the stan- 
dard model. For example: Are the TLSs uniformly dis- 
tributed in space, or do they preferentially appear at 
boundaries of clusters of atoms/molecules, as suggested 
in [l3l | . Another open question is the nature of the inter- 



action between the TLSs and the SM e.g., is it dipolar 
as suggested in 0? A method to obtain this important 
information using SM spectroscopy was suggested in Q 
(see details below). 

A theoretical investigation of the distribution of ran- 
dom line shapes of SMs in glass was carried out in p^[l5| 
based on the standard model approach . Interestingly, 
the theoretical results obtained in |3 showed that the 
problem of random line shapes of SMs in glasses is related 
to Levy statistics, thus the generalized Central Limit 
Theorem applies to this problem. This connection 
to Levy statistics is a manifestation of long-range inter- 
actions between the TLSs in the glass and the SM (see 
details below). We note that Levy statistics is known to 
describe several other long-range interaction models in 
diverse fields such as turbulence and random mag- 
netic systems Stoneham's theory of inhomoge- 
ncous line broadening in crystals with defects is based on 
long-range forces, and it can be interpreted in terms of 
Levy statistics |15| . 

In this paper we analyze the statistical properties of 
random line shapes in aglass and compare it to the the- 
oretical predictions in [lj|- In Fig- Qwe show eight lines 
of single tetra-tert-butylterrylene chromophores embed- 
ded in an amorphous polyisobutylene matrix at T = 2 K. 
For experimental details see [23. The lines are typically 
multi-peaked similar to the numerical predictions in 9|. 
Each line is different from any other line, since each indi- 
vidual SM is in a unique environment. The multi-peaked 
behavior of the line shapes in Fig. ^can be qualitatively 
explained using standard model arguments. If a SM is 
coupled to a single slow flipping TLS embedded in its 
vicinity, we expect that when the TLS flips from its up 
state to its down state or vice versa, the SM absorption 
frequency will shift. In this case the line shape of the SM 
is a doublet. It follows that the frequency of a SM cou- 
pled to N such independent two-level systems will jump 
between 2^ states. Hence as shown in Fig. ^ the single 
molecule's line shape will be composed of 2^ peaks. The 
width of these peaks depends on the dynamics of TLSs 
situated far from the molecule, fast dephasing processes. 
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and the life time of the electronic transition |9l lT3.l2lLl2^ . 

To obtain the lines in Fig. Qwe used the spectral trail 
technique introduced by Moerner and co-workers [23l |. 
following the spectral activity of the molecule during a 
scan time which in our experiments was fixed to be 120 
Sec. Following the jump history of the SM enables us to 
identify the peaks in Fig. ^ originating from a SM. 
Without using this method it is practically impossible to 
say if the lines in Fig. ^ are due to contribution from 
several SMs or originate from a SM. 
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FIG. 1: Line shapes of single tetra-tert-butylterrylene chro- 
mophores embedded in an amorphous polyisobutylene matrix 
at T = 2 K. Note the doublet and quartet features of some 
of the line shapes due to strong coupling to one or two TLSs. 
One of the molecules has three peaks, indicating the possibil- 
ity that this molecule is coupled to a three-level system. This 
rare type of behavior is not consistent with the standard tun- 
neling model. Another possibility is that measurement time 
was not long enough in this case. 

Random multi-peaked lines as shown in Fig. ^ ^^'^ ^ 
novel feature of SM spectroscopy, not observed in en- 
semble averaging techniques. To obtain information on 
the SM-glass system we investigate the distribution of 
the line shapes of SMs in a glass. Mathematically we 
are dealing with the question of the distribution func- 
tion of a function [i.e., the line J(a; r , ) is a function of 
laser frequency]. Recently, Refs. suggested to 

characterize the line shape of each molecule by its cu- 
mulants ki, K2, K3, .... Thus each line is characterized by 
an infinite set of cumulants, which are random variables 



that vary from one molecule to the other, thus reflect- 
ing the disordered nature of the glass. The distribution 
functions P{ki), P{k2), P{k3), ■ ■ ■ completely character- 
izes the statistical properties of the line shapes of SMs in 
a glass. The cumulants are obtained from the moments 
of the line shape 24], m„ = J lu^I {lol) dwi, according to 
the well-known relations ki = rni, K2 = m2 — {niiY etc. 
From 244 SM spectra, we obtained the histograms of the 
first two cumulants, P(ki), P{k2), which yield important 
information on the glass-SM system. We consider here 
the distribution of the cumulants and not the distribu- 
tion of the moments, since the cumulants were predicted 
to be described by Levy statistics |l4j . 

To better understand the meaning of our data analy- 
sis and the relationship to Levy statistics, it is useful to 
recall five main assumptions of the model used in 0, 0] . 
(i) The absorption frequency of the SM follows a stochas- 
tic trajectory described by u}{t) = + J2n=i iri{t)i'n{r), 
where wq is the bare frequency of the molecule. ^n(t) 
are random functions of time which follow a two-state 
process, ^n(^) = 1 when the n-th TLS is in state up or 
^n{t) — when it is in state down. Thus the fiips of the 
TLSs induce spectral diffusion. The flipping rate between 
the up and down states is determined by a rate R which 
varies from TLS to TLS [the distribution of jumping rates 
of TLSs in a glass spans many orders of magnitude from 
nano seconds to (at least) days], (ii) The SM frequency 
shifts are: f„ = 27ra^ (il) ^(^rjg-- The most important 
ingredient of the theory is the long-range interaction be- 
tween the TLS and SM, v cx reflecting the assump- 
tion of dipolar interaction. This long-range type of inter- 
action is the first important ingredient in the relationship 
between SM spectroscopy in glasses and Levy statistics. 
Other parameters controlling the frequency shifts are: a 
the SM-TLS couphng constant, E = \/A^ + ,P the en- 
ergy splitting of the TLS, ^ (il) a dimensionless func- 
tion of order unity describing the orientation of the TLS 
and SM, and finally the random parameters of the TLS: 
A (asymmetry parameter) and J (tunneling matrix ele- 
ment), (iii) The TLSs are uniformly distributed is space 
and are non-interacting. This assumption is the second 
important condition for Levy statistics to be valid, (iv) 
The standard tunneling model is valid; this model deter- 
mines the distribution of parameters describing Vn and 
Rn as well as the density of TLSs. Note, however, that 
the Levy statistics and results in Ti^ are not limited to 
this model, (v) The stochastic Kubo-Anderson theory of 
line shapes is applicable, implying weak laser fields. Un- 
der these conditions the following two limiting behaviors 
were found ^14] . 

The first case corresponds to the fast-modulation limit 
Vn ^ Rn for all TLSs in the vicinity of the molecule. In 
this case, also called motional narrowing limit, the lines 
of individual molecules are Lorentzian in shape. Then the 
lines are characterized by two parameters only, e.g., the 
width at half maximum and the center location. The dis- 
tributions of these two parameters are Levy stable laws 
[1^ . From Fig. ^ it is clear that the fast modulation 
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limit does not describe our experimental results. 

The second case corresponds to the slow-modulation 
limit. If Rn ^ i^n for all the TLSs in the vicinity of the 
molecule, the shape of the line is random and typically 
multi-peaked . In this slow-modulation limit the dis- 
tributions of line shape cumulants, P{ki), P{k2), P{k3), 
etc, are all Levy stable. Specifically, the probability den- 
sity function of the first cumulant ki is given by the 
symmetrical Levy density, P{ki) = ^loC^^i); namely the 
Lorentzian 
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where zi is a scaling parameter which can be calculated 
from the theory in . In Fig. El we show that our 
experimental results are compatible with the theoreti- 
cal prediction. We also fitted our results to a Gaussian 
probability density (not shown) and found that the dis- 
tribution of the first cumulant is definitely not Gaussian. 

We note that the reference frequency, determining the 
laser detuning in Fig. ^ was choosen on the maximum 
intensity of the spectrum of the SM ujmax- While the 
theoretical reference frequency in |l4j was the bare fre- 
quency of the molecule ujn. Using numerical simulation 
based on the approach in [9j , we observed that also when 
the reference frequency is choosen as tOmax distribution of 
first cumulant is well fitted by a Lorentzian, in agreement 
with the experimental results in Fig. 12 The value of zi in 
the two approaches slightly differs, as we discuss below. 
We now discuss the distribution of the second cumulant 
(i.e. the variance) which is not sensitive to definition of 
the reference frequency. 

The distribution of the second cumulant K2 is given 
by the one-sided Levy stable law, P{n2) — ^i/2,i('*2), 
namely Smirnov's one-sided probability density 



P{n2) = 




-3/2 



exp 



^1/2 \ 



(2) 



where the scaling parameter zi/2 was derived in [T3 |. As 
shown in Fig. |3 the experimental data for the distribu- 
tion of K2 are compatible with the theoretical prediction; 
the long tail of the Levy stable law is visible. Yet, data 
of a larger number of molecules is needed to improve the 
statistical fiuctuations. The analysis of a larger number 
of molecules will also enable us to compare theory and 
experiments for the higher-order cumulants K3, K4 etc. 

As mentioned, the Levy behavior is due to the long- 
range dipolar interaction of a SM with many TLSs, 
hence the Gauss or Levy Central Limit Theorem ar- 
guments are expected to hold. The Levy central limit 
theorem applies since averaged frequency shifts diverge; 
(i^) cx /g r"^^^ /r^dr = 00 where d = 3 is the dimension- 
ality of the problem (for mathematical details see |0|). 
Hence, the Levy behavior obtained for the distributions 
P{ki) and P{k2) is used to test the assumptions of dipo- 
lar interactions and uniform-distribution of the TLSs in 
space. The information about the random distribution of 



the parameters of the glass, i.e. A and J, are contained 
in the values of Zi/2 and zi which we now discuss. 




FIG. 2: Probability density of the first cumulant ki (units 
are GHz). The dots are experimental results, the curve is a 
one-parameter fit to a Lorentzian. 

There are two distinct types of parameters describing 
the SM-glass system: Those describing the bath of TLSs 
and the coupling constant a which depends on the prop- 
erties of the SM probe. An important unsolved prob- 
lem is how to extract information from SM experiments 
which is sensitive only to the statistical properties of the 
TLSs and is not affected by a The scaling param- 
eters Zl and zi/2 depend on a, the magnitude of these 
two parameters depends also on the precise modeling of 
the orientation function ^'(O) Thus zi and 2^/2 are 
not universal functions, in the sense that they depend 
on properties of the SM under investigation and not on 
the properties of the glassy state (which are supposed to 
be universal, according to standard tunneling model). In 
fact, Donley et al ^ suggested that the coupling constant 
a itself should be a random variable. This may seem to 
limit our ability to investigate low-temperature glasses 
with SM spectroscopy. If the statistical analysis of line 
shapes depends on an additional unknown distribution 
function of the coupling constant (besides the standard 
distributions of the glass parameters), fitting of data to 
the theory becomes rather arbitrary. 

However, based on Eq. 9 in 14] one can show that the 
ratio 2:1/2/2:1 depends only on the statistical properties 
of the glass and not on the distribution of the coupling 
constant a. More precisely. 
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where the averaging is performed over the TLS parame- 
ters A and J [iJI- Since Eq. Q is independent of the 
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exact distribution of a, it is a useful tool for describing 
the behavior of glasses. To derive Eq. Q we assumed 
that the random variable a is independent of the glass 
parameters A and J. 
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FIG. 3: Probability density of the second cumulant K2 (units 
GHz^). The dots are experimental results, the curve is a one- 
parameter fit to Smirnov's probability density. 

From our fits we find Zi/2 — 0.175 GHz and zi = 
0.0485 GHz, which yields zi/2/21 = 3.6. The theoret- 
ical prediction, based on Eq. Q yields Z112/ zi = 2.4. 
The deviation between theory and experiment is now 
explained. As mentioned, the theoretical prediction is 
based on the assumption that the bare frequency of the 
molecule wq is the reference frequency for measurement. 
Using numerical simulations with the parameter set 
relevant for our experiment (details will be published) 
we find that the ratio of zi when loq is the reference fre- 
quency and zi when Wmax is the reference frequency, is 
z"" /z^"""^ — 1.6. Varying the value of a in our simula- 
tions, in the interval lOGHznm^ < a < 40GHznm'^ we 
observed that the ratio z^'Y^^j^""^ does not depend on 



a. As expected 2;"° > 2;""""= since the value of ki be- 
comes smaller (in statistical sense) if we assign the origin 
to the maximum of the spectrum. Using the correction 
factor zp/z"""""^ = 1.6 we find that the theory yields 
•2^1/2/^1 = 3.8. Taking into account that the standard 
tunneling model does not address the chemical compo- 
sition of the disordered system or the chemical and geo- 
metric details of the SM under investigation, we believe 
that the theoretical result is in surprisingly reasonable 
agreement with experiment. Measurements of the ratio 
zxjijzx for other types of SMs and glasses and for wider 
range of control parameters (i.e., temperature and scan 
time) will show whether SM data are compatible with 
the universal predictions of the standard model. 

Our approach based on the analysis of the distribu- 
tions of line shape cumulants is different from a second 
approach used in other studies 0, |H ■ Previously, the 
line was characterized by its width at half maximum and 
the distribution of the random line shapes was described 
by distribution of line widths. Clearly such a method 
does not capture the multi-peaked behavior of the SM 
lines shown in Fig. ^ An approach based on the dis- 
tribution of line widths is useful in the fast-modulation 
limit, which as mentioned, (i) obeys Levy statistics but 
(ii) is not relevant to the experiment under consideration. 

To summarize SM spectroscopy is an excellent method 
to investigate disordered systems by removing the ensem- 
ble averaging found in other techniques. We have shown 
that our experimental results are compatible with Levy 
statistics and with standard tunneling model predictions. 
In particular the following two assumptions are reason- 
able: (i) The two-level systems are uniformly distributed 
in space, (ii) The frequency shifts are caused by dipolar 
interactions v ^ We introduced the universal ra- 
tio, Z1/2/ zi which is sensitive to details of the standard 
model, but not to the coupling of the SM to the TLSs 
in the glass (i.e., not to a). The comparison between 
the theoretical and the experimental value of this ratio 
can be used to test the validity of the standard model 
predictions. 
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